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Abstract A qualitative analysis is given to the data on the full magnetic and 
velocity vector fields in a growing sunspot group, recorded nearly simultaneously 
with the Solar Optical Telescope on the Hinode satellite. Observations of a young 
bipolar subregion developing within AR 11313 were carried out on 9-10 October 
2011. Our aim was to form an idea about the consistency of the observed pattern 
with the well-known rising-tube model of the formation of bipolar active regions 
and sunspot groups. We find from our magnetograms that the distributions of 
the vertical [i?v] and the horizontal [Rh] component of the magnetic field over 
the area of the magnetic subregion are spatially well correlated; in contrast, the 
rise of a flux-tube loop would result in a qualitatively different pattern, with the 
maxima of the two magnetic-field components spatially separated: the vertical 
field would be the strongest where either spot emerges, while the maximum 
horizontal-field strengths would be reached in between them. A specific feature, 
which we call the bordering effect, is revealed: some local extrema of are 
bordered with areas of locally enhanced This effect suggests a fountainlike 
spatial structure of the magnetic field near the By extrema, which is also hardly 
compatible with the emergence of a flux-tube loop. The vertical-velocity field in 
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the area of the developing active subregion does not exhibit any upflow on the 
scale of the whole subregion, which should be related to the rising-tube process. 
Thus, our observational data can hardly be interpreted in the framework of the 
rising-tube model. 

Keywords: Velocity Fields, Photosphere; Magnetic fields; Active regions; Sunspots| 


1. Introduction 


In the developing active regions (ARs), subphotospheric magnetic helds and fluid 
motions are strongly coupled according to the laws of magnetohydrodynamics. 
Obviously, the new magnetic flux emerges to the solar surface from deeper layers; 
if the emerging field is sufficiently strong, sunspot groups may develop. However, 
an intriguing issue related to this process is whether the strong magnetic field 
plays an active, primary role, having been prepared in the depth of the convection 
zone and then emerging through the visible surface of the Sun, or the primary 
role is played by the fluid motion, which amplifies and structures the initially 
moderate or weak magnetic field. 

The first of these alternatives is typically represented by the popular rising- 
tube model (RTM) according to which the magnetic field of a bipolar sunspot 
group originates from the emergence of an fl-shaped loop of a coherent flux tube 
of strong magnetic field. As usually assumed, the flux tube forms in the general 
toroidal^ solar magnetic field deep in the convection zone, and the field that this 
tube carries upward is already strong before the r ise. T his mechanism received 
much attention after a well-known study by Parker (119551 ). who invoked magnetic 
buoyancy to account for the rise of flux-tube loops. Later, the RTM has been 
revisited by a number of investigators over several decades. Interesting consider¬ 
ations of this model were suggested, in particular, by Caligari, More no-Ins ertis, 
and Schiissler (|l995|) and Caligari, Schiissler, and Moreno-Insertis (|l998l l: nu¬ 
merical simulations of this mechanism based on full systems of M HD e quations 
have also bee n carr ied out (see, e.g., Fan, Featherstone, and Fang. l2013l : Rempel 
and Cheung. 12014 . and references therein). 

The concept of RTM, which agrees with such important regularities of solar 
activity as Hale’s polarity law and Spbrer’s law of sunspot-formation latitudes, 
appeared to be very attractive in the epoch of moderate capacities of obser¬ 
vational instrumentation. For this reason, the idea of RTM remained virtually 
indisputable for a long time. Currently, the RTM still receives attention in both 
analyses of observational data and numerical simulations. In particular, under 
the as sumpt ion that a twisted fl ux tu be (rope) rises in a developing AR, Luoni 
et al. ( 201ll l and Poisson et al. ( 2015ll estimate the magnetic helicity using the 
configurations of magnetic polarities in observed ARs. 

However, some implications of the RTM were found to be hardly compatible 
with observations. In particular, if a flux-tube loop is rising, two very impressive 


^As frequently done in the literature on stellar and planetary dynamos, we use here the terms 
toroidal and azimuthal as synonyms, although they are not mathematically equivalent. 
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manifestations of this process should definitely be observed but, as we shall see, 
they do not conform with the observations discussed below. 


i) As the loop is rising, strong horizontal magnetic fields should emerge on the 
scale of the entire AR. It is the pattern of magnetic fields in a developing 
AR that will be the subject of our discussion here, and we shall find no 
manifestations of the loop rise. 

ii) Intense spreading from the site of the loop emergence should be observed 
over th e who le AR. This is actually not the case; in particular, Pevtsov and 
Lamb ( 200fih “observed no consistent plasma flows at t he fut ure location of 
an active region before its emergence” and Kosovichev ( 2009l l finds “no evi¬ 
dence for large-scale flows indicating future appearance a large-scale magnetic 
structure without signs of a large-scale horizontal flow”. A further example 
of the horizontal-velocity pattern without spreading on th e scale of the whole 
developing AR is given by Getling, Ishikawa, and Buchnev ( 20I5ll (hereinafter. 
Paper I). 


In addition, it is worth mentioning two other important considerations that are 
also not in favour of the RTM: 


iii) If the RTM is adopted, one has to account for the origin of the coherent 
tube of strong magnetic field deep in the convection zone. Various assump¬ 
tions have been made to this end, which differ in their plausibili ty and the 
appropriateness of their starting points (see, e.g., a review by Fan, 20091 . and 
references therein). Furthermore, before the emergence on the photospheric 
surface, such an intense flux tube should affect the structure of the convective 
velocity field, which also is not actually observed. 

iv) Joy’s law of the latitudinal dependence of the tilt angle of bipolar sunspot 
groups is not sensitive to the amount of the emerging magne tic flu x in contrast 
to what coul d be expected: as Kosovichev and Stenflo ( 2008[l and Koso¬ 
vichev ( 2009ll write, their “new statistical study of the variations of the 
tilt angle of bipolar magnetic regions during the flux emergence questions 
the current paradigm that the magnetic flux emerging on the solar surface 
represents large-scale magnetic flux ropes (fl-loops) rising from the bottom 
of the convection zone” (italicised by us. - Authors). 


These are the most important points of doubt about the universal applica¬ 
bility of the RTM; we shall not go here into further details, since we already 
discussed some of them in Paper I. In view of the contradictions of the RTM 
with observations and difficulties of accounting for the origin of the intense flux 
tube and some other features of the process, the rising-tube mechanism no longer 
appears to determine a paradigm in the studies of the development of ARs. 

As alternatives to the RTM, various mechanisms of in situ magnetic-field am¬ 
plification and structuring have been suggested. Among them, MHD mechanisms 
of inductive excitation of magnetic fields strongly coupled with fluid motions are 
usually referred to as local dynamos. This term is frequently used only in the 
context of small, granular scales and under; in contrast, we associate the concept 
of local dynamo with a wider range of scales, including the sizes of the whole 
ARs (mesoscales). 
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The idea of local MHD dynamo traces back to Gurevich and Lebedinsky 


(|l946f l. who related the amplification process to the effects of plasma motions; 
however, they did not attribute these motions to convection and even did not 
specify any particular type of motion. Tverskoi ( IQfifih used a very simple model 
to demonstrate that the magnetic field can locally be amplified and structured 
by cellular magnetoconvection. Later, numerical simulations revealed the role 
of local dynamos as the producers of fairly disordered, intermittent magnetic 
fields on v ery sm all (down from t he gra nular) scales ('Cattaneo. ll999l : Vdgler and 
Schiissler, 2007 : Kitiashvili et ai, I 2 OI 5 I. etc.) In essence, only Stein and Nordlund 
(|2012n used “realistic” numerical simulations to describe a convective mechanism 
capable of producing mesoscale amplified magnetic fields. They investigated the 
formation of an AR via the flux rise due to convective motions in the upper 
portion of the convection zone. The computed scenario does not imply the pre¬ 
existence of a coherent flux tube. A uniform, untwisted, horizontal magnetic field 
is initially present, and magnetic loops subsequently form over a wide range of 
scales. However, the initial field is required to be relatively strong, and only 
a moderate magnetic-field amplification (by a factor of about three) can be 
achieved. Previously, we discussed local convective dynamos in more detail (see 
Paper I). 

Sunspot-formation mechanisms differing from local dynamo were also sug¬ 
gested. Kitchatinov and Mazur ( 200Cll l investigated a hydromagnetic instability 
that can act on scales large compared to the granular size, producing a magnetic- 
flux concentration similar to those observed in sunspots. The process crucially 
depends on the presence of fluid motion and on the quenching of eddy diffusivity 
by the enhanced magnetic field with the plasma cooling down. This is a local 
mechanism, which, however, is not a dynamo in the strict meaning of this term. 

Another local mechanism, which does not qualify as a dynamo, is related 
to the so-called negative-effective-magnetic-pressure instability (NEMPI), see 
Warnecke et al. ( 2013 . 2015^ and references therein. It results from the suppres¬ 
sion of the total turbulent pressure (the sum of hydrodynamic and magnetic 
components) by the magnetic field. 

We present here a qualitative consideration of some observational data from 
the standpoint of items (i) and (ii) in the above list to decide whether or not these 
data can be interpreted in terms of the RTM. To obtain relevant data, we devel¬ 
oped an observational program of studying the evolution of both the velocity and 
magnetic fields in growing ARs. This program (operation plan) is intended for 
implementation with the Solar Optical Telescope (SOT) on the Hinode spacecraft 
(Tsuneta et al., 20081 : Suematsu et ai, 20081 : Shimizu et ai. I 2 OO 8 I and has been 
designated as HOPI8I (http://www.isas.jaxa.jp/home/solar/hinode_op/hop.php? 
hop=0181). It consists in simultaneously recording and analysing the dynamics of 
the full-vector velocity and magnetic fields on the photospheric level. Previously, 
we presented some preliminary results of our study (Paper I). New features of 
the AR evolutionary pattern inferred from our data will be given here. 
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2. Observations and Data Processing 

A bipolar magnetic structure, which emerged within AR 11313, was observed 
at its early evolutionary stage, on 9-10 October 2011; the AR was then near 
the centre of the solar disc. Five 2-h-long observational sessions were carried out 
with intervals that varied from 3 h 40 min to 6 h 30 min (see a summary in 
Table HI). 


Table 1. Summary of observational sessions 


Session 

No. 

Date 

Session 

meantime 

Interval 
from (1) 

1 

09 Oct. 2011 

19:31:15 

00:00 

2 

10 Oct. 2011 

01:06:16 

05:35 

3 


06:30:15 

10:59 

4 


15:00:14 

19:29 

5 


21:46:15 

26:15 


During each session, a 150" x 163" field of view (FOV) was observed using 
the Narrowband Filter Imager (NFl) of the SOT at two wavelength positions of 
Fel A 5776 A with a time cadence of 2 min and a pixel size of 0.16". This yielded 
a series of photospheric images, which can be used to calculate horizontal ve¬ 
locities [uh], and a series of Dopplergrams representing the line-of-sight (LOS), 
or vertical, velocities [uv]- Simultaneously, t he sam e FOV was scann ed with the 
Spectra-Polarimeter (SP; see Ichimoto et al. 12008 : Lites et al., 2013h one or two 
times a session. The SP scan was done in the so-called fast mode with a pixel 
size of 0.32", taking 32 min to obtain one SP map. To derive full-vector magnetic 
fields from these SP observations, we used the MERLIN inversion code (Lites 
et al., 2 OO 7 II . which assumed a Milne-Eddington atmosphere. 

The processing of the data included: 


i) Subsonic filtering based on Fast Fourier Transform; 

ii) Constructing Dopplergrams; 

iii) An intensity-scaling procedure enhancing the image contrast by means of cut¬ 
ting off the tails of a pixel-intensity histogram and subsequent linear mapping 
of the remaining portion of the histogram onto the whole admissible intensity 
range; 

iv) Alignment of the magnetograms corrected for irregularities of the SP scanning 
process, with the properly rescaled images and Dopplergrams obtained at 
the meantime of the scan (the correction was done by choosing 30 or more 
reference points, typically pores or easily identifiable fine details of spots, in 


^ Since the area of interest was located near the solar-disc centre and, moreover, corrections 
for projection effects are not important from the standpoint of our goal, we do not make 
difference here between the line-of-sight and the vertical component and also between the 
transversal (tangential) and horizontal vector components. 
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both the FG image and the image accompanying the magnetogram, and sub¬ 
sequently bringing them into coincidence by means of affine transformations); 

v) Determination of the horizontal-velocity field using a technique based on the 
same principle as the standard method of loc al cor relation tracking (LCT) 
but more reliable (see Getling and Buchnev, 201Cll . for a description) and 
construction of cork-trajectory maps. Our technique differs from the standard 
LGT procedure in a special choice of trial areas (“targets”), whose displace¬ 
ments are determined by maximizing the correlation between the original 
and various shifted positions of the target. Specifically, an area is chosen as a 
target in a certain neighborhood of each node of a predefined grid if either the 
contrast or the entropy of the brightness distribution reaches its maximum 
in this area. The horizontal velocities obtained are then interpolated to the 
positions of imaginary “corks” using the Delaunay triangulation and afhne 
transformations specified by the deformation of the obtained triangles at the 
time step considered; 

vi) Elimination of the Sun’s rotation from the fields of LOS velocities; 

vii) Reducing the mean LOS velocity to zero in each map. 


The MERLIN code yields the magnitude, inclination and azimuth of the 
magnetic-field vector (in parallel with a continuum image), which we then con¬ 
vert into the LOS and tangential components. In addition to the original maps 
of the vertical velocity and magnetic-field components, we used smoothed maps 
of these quantities. 


3. Results 

3.1. Evolution in White Light 

The photospheric images of AR 11313 obtained during the five observational ses¬ 
sions are shown in Figure [T] They cover a time interval somewhat exceeding one 
day (see Tabled]). It can be seen that a bipolar sunspot group has already formed 
by the first observational session. At nearly the same time, a new, minor sunspot 
group starts developing between the main spots of the previously developed 
group, in the left half of the FOV. This process becomes mainly accomplished 
by the third session, after which the spots of the newly born group continue 
growing and its leading and trailing spots move apart. 

3.2. Magnetic-Field Patterns. Bordering Effect 

We consider here some features of the developing magnetic and velocity fields 
deducible from Figures [SHU which refer to sessions (l)-(5), respectively. The 
left-hand panel of each figure makes it possible to compare between the origi¬ 
nal (detailed) maps of the vertical [B^] and horizontal [i?h] components of the 
magnetic field. As for confronting the vertical magnetic and velocity helds, a 
specific noise of fine details in the original maps makes comparisons between 
them difficult. In view of this, we use smoothed magnetic and velocity fields for 
such comparisons (right-hand panels). 
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Figure 1. Evolution of AR 11313 in visible light. Images (a)-(e) were taken at the mean times 
of observational sessions (l)-(5) (see Table [Til. 
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Figure 2. Left: comparison between the vertical (grey-scale map) and horizontal (con¬ 
tours) components of the magnetic field for observational session (1); the Bv range is 
[—1941 G, 2477 G]; contour levels for are 800 G (red) and 1400 G (yellow). Right: com¬ 
parison between the patterns of the smoothed vertical velocity field (grey-scale filled contour 
map with the dotted contours for zero velocity) and smoothed vertical magnetic field (colour 
contours) for the same session; the contour increment is 800 G for magnetic field and 0.2 km s“^ 
for velocity; red contours correspond to Bv > 0, green dot-dashed contours to Byj = 0 and 
yellow contours to B^ < 0; the size of the smoothing window is 3.5 Mm for magnetic field and 
7.4 Mm for velocity. The dark areas in both grey-scale maps correspond to negative values 
(vectors directed upward) and light areas to positive values (vectors directed downward). 




Figure 3. Same as in Figure[2]but for session (2); the Sv range is [—2121 G, 2460 Gj. 


Let US first note a remarkable feature of the magnetic field, which we call the 
bordering effect. It can be revealed in the magnetograms of all the five sessions 
by examining the left-hand panels but is most pronounced in Figures [S] and [51 
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Figure 4. Same as in Figure[2]but for session (3); the Bv range is [—2112 G, 2467 G]. 
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Figure 5. Same as in Figure[2]but for session (4); the range is [—2177 G, 2217 Gj. 




While most local (mainly positive) extrema of spatially coincide with areas of 
low i3h, these areas are partially bordered with arched areas of locally enhanced 
Sh. 

The bordering effect admits a qnite straightforward interpretation. Assume 
that magnetic field lines form a bundle issuing from the subphotospheric layers 
and diverge in its top section like the water jets in a fountain. In this case, 
will obviously reach its peak magnitudes in the central part of the bundle, while 
Bh will be maximum in an annulus surrounding the bundle. 

A common feature of all the magnetograms is fairly good, visually noticeable 
spatial correlation between the distributions of By and B^ (see left-hand panels 
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in Figures [5Hn]). However, for a given sign of Hv, its pattern typically exhibits a 
relatively uniform spatial shift with respect to the corresponding, similar pattern 
of i?h (the contour representation of clearly visualises this shift only in the 
areas of < 0). The shifts for the two signs may differ in their magnitude and 
direction; they do not exceed 2 Mm. 

The shift of the i?h pattern relative to the B^ pattern may be indicative of a 
systematic increase in the tilt of the magnetic field lines in the direction of this 
shift. This feature thus proves to be akin, in a sense, to the bordering effect, 
being a manifestation of the three-dimensional structure of the magnetic field. 

To make further inferences from the magnetograms, we have to remember 
a feature that should be expected in the case of the flux-tube rise and was 
mentioned above as item (i) in the list of RTM doubtful points (see Introduc¬ 
tion) . Specifically, strong horizontal fields should connect the emergence areas of 
the main spots, forming elongated features where the horizontal-field strengths 
values are comparable with the vertical-field strengths observed in the spots. 

As for our magnetograms, they display signatures of finely structured hor¬ 
izontal field in the form of very narrow elongated features stretching between 
the main magnetic poles. These features are present in the obtained i?h distri¬ 
butions for sessions 3, 4 and (in a very faint appearance) 5 but can hardly be 
distinguished in their contour representations in the respective figures. In this 
context, the following should be noted. First, these features became noticeable 
after the formation of the bipolar configuration of strong magnetic field but 
were completely absent at the stage of development of this configuration. In 
contrast, if a flux-tube loop rose, we would see stretched features (not necessarily 
finely structured) earlier, during the sunspot-group development. Second, as the 
magnetograms demonstrate, the horizontal field in these features is considerably 
weaker (a factor of two to three) than the peak Hh values, which, in turn, do not 
exceed the peak values of B^ (see Tabled). Therefore, the scenario of magnetic- 
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field development is not consistent with the RTM-based expectation indicated 
as item (i) in the list of RTM doubtful points (see Introduction). 


Table 2. Characteristic values of Bv and Bh in the developing 
AR 


Session 

No. 

Bv range 

[G] 

B\^ range 

[G] 

Typical values 

in fine features [G] 

1 

(-1900)-2500 

0-1900 

800-1200 

2 

(-2100)-2500 

0-1800 

600-1000 

3 

(-2100)-2500 

0-2100 

700-1100 

4 

(-2200)-2200 

0-1800 

700-1100 

5 

(-2400)-2300 

0-2200 

600-1100 


3.3. The Magnetic vs. the Velocity Field 

Now let us compare the smoothed and Uv fields shown in the right-hand 
panels of Figures and discuss the relationship between them. It can be seen 
from our velocity maps that no marked upflow can definitely be related to the 
entire area where the rise of a flux-tube loop should be expected according to the 
RTM. In Figured! several moderate, localised upflows are present in this area, 
while the strongest upflows are present not far from the lower left corner of the 
FOV and near the existing sunspots. In Figured! an upflow-velocity maximum 
is located below the area of interest, being comparable in its magnitude with 
several other local maxima. Nothing similar to a strong upflow associated with 
the hypothetical rising tube can be seen in Figured! not to say about Figures [5! 
and in! 

It is also worth remembering that the Uh held that we constructed based 
on the same observations as those discussed here (see Paper I) exhibits neither 
any spreading how on the scale of the whole developing group nor any hows 
qualitatively different from normal mesogranular and supergranular convection 
unaffected by any large-scale disturbances. 

Thus, the observed patterns of both the vertical and the horizontal velocity 
in the developing active subregion are at variance with the RTM-based expec¬ 
tation mentioned above as item (ii) in the list of RTM doubtful points (see 
Introduction). 

It is also interesting to note the following particular feature of the Uv held in 
the vicinity of the leading (rightmost) spot of the large group that was formed 
before the smaller group started developing. If we compare Figure [T^-c with 
and Uv distributions in the right-hand panels of Figures [5H3! we can hnd that, on 
two sides of the leading spot, there are an uphow (above the spot in the FOV) 
and a downhow area (below it). This feature can also be distinguished (although 
is much less pronounced) in Figures [5! and [5! 
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4. Discussion and Conclusion 


Our analysis of the observational data for the development of the minor sub- 
region within AR 11313 on 9-10 October 2011 reveals a noticeable discrepancy 
between the observed development pattern of the magnetic and velocity fields 
and the RTM-based expectations. 

First, we have found that the distributions of and Rh over the area of the 
growing magnetic subregion are spatially well correlated, with a shift between 
the entire patterns of and B^ by a distance of no larger than about 2 Mm. 
The rise a flux-tube loop should result in a qualitatively different pattern. The 
maxima of the two magnetic-field components would be spatially separated in 
this case: the vertical field would be the strongest where either of the main 
spots emerges, with the maximum horizontal-field strengths reached in between 
them. Moreover, the horizontal field would form elongated features connecting 
the main spots, being comparable with the vertical-field strengths in these spots. 
In contrast, we observe finely structured elongated features, which emerged no 
earlier than the bipolar configuration of strong magnetic field had completely 
formed and were considerably weaker than the vertical field in the spots. 

The feature that we call the bordering effect provides additional evidence 
against the rise of a tube, since it demonstrates a fountainlike three-dimensional 
structure of individual local magnetic-field maxima. It becomes pronounced by 
the final formation stage of the active subregion under study (Figure 0]). Such a 
structure is hardly compatible with the emergence of a whole tube loop, which 
should produce local maxima of the vertical magnetic field only at the feet of 
the loop. 

Second, the flow pattern observed in the area of the developing active sub- 
region is not consistent with the idea of the flux-tube-loop emergence. There 
is no upflow on the scale of the whole subregion, which should be related to 
the rising-tube process. In addition, the horizontal-velocity field (considered in 
Paper I) does not exhibit any spreading flow on the scale of the entire growing 
magnetic region; instead, some motions resembling normal mesogranular and 
supergranular flows appear to b e pres e rved. In this context, it is worth noting 


that, many years ago, Bumba ( 1963 . 1967ll and Bumba and Howard ( 1965 1 


found that the growing magnetic fields do not break down the pre-existing 
convective-velocity field but come from below “seeping” through the network 
of convection cells. The development of a sunspot group seems to be controlled 
by the supergranular network, and the lines of force of the strong local magnetic 
fields are nearly collinear with the streamlines of the photospheric plasma. This 
can naturally be understood if the magnetic field is assumed to be formed by 
convective motions. 

Thus, our observational data can hardly be interpreted in the framework of the 
rising-tube model. However, our inferences are based on a single AR-emergence 
event, which is insufficient to assess how typical the observ ed sce nario is. 

The emergence of 41 AR is analysed by Poisson et al. ( 20151 ) based on ob¬ 
servations of the vertical magnetic fields. They use a model of the twisted flux 
tube (rope) to describe the time variations in the twist (and, accordingly, in the 
magnetic helicity). However, as can be observed, the assumption of the flux-tube 
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rise is the starting point of their investigation and its results are discussed in 
comparisons with numerical simulations of the rising tubes. Since no possibilities 
lying beyond the framework of this concept are considered, it is not obvious that 
such alternatives should be ruled out in the cases considered. 

The issue of the comparative role of different really possible active-region gen¬ 
eration mechanisms thus remains open, and further observations are needed. In 
our view, they would be most promising if they include recording both magnetic 
and velocity fields. We plan such observations for the near future. Since the 
rising-tube mechanism does not appear to be universal, local convective dynamo 
as the producer of the magnetic fields of active regions deserves close attention 
(see Paper I for a discussion of the relevance of this mechanism). 
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